


207150 


Surface Turbulent Fluxes Over Pack Ice 
Inferred from TO VS Observations 

R. W. Lindsay*, J. A. Francis**, P. O. G. Persson***, 

D. A. Rothrock*, and A. J. Schweiger* 

* Polar Science Center 
Applied Physics Laboratory 
University of Washington 
1013 NE 40th Street 
Seattle, WA 98105 USA 



//j - ' V- 

<33 7 2.-^7 


** Institute of Marine and Coastal Sciences 

P.O. Box 231, Rutgers University 

New Brunswick, NJ 08903-0231 USA 

*** Cooperative Institute for Research in Environmental Sciences 
University of Colorado / NOAA, Environmental Technology Laboratory 
325 Broadway 
Boulder, CO 80303 USA 


Abstract 

A one-dimensional, atmospheric boundary layer model is coupled to a ther- 
modynamic ice model to estimate the surface turbulent fluxes over thick sea 
ice. The principal forcing parameters in this time-dependent model are the air 
temperature, humidity, and wind speed at a specified level (either at 2 m or at 
850 mb) and the downwelling surface radiative fluxes. The free parameters are 
the air temperature, humidity, and wind speed profiles below the specified 
level, the surface skin temperature, the ice temperature profile, and the surface 
turbulent fluxes. The goal is to determine how well we can estimate the turbu- 
lent surface heat and momentum fluxes using forcing parameters from atmo- 
spheric temperatures and radiative fluxes retrieved from the T1ROS-N 
Operational Vertical Sounder (TOVS) data. 

Meteorological observations from the Lead Experiment (Lead E x, April 
1992) ice camp are used to validate turbulent fluxes computed with the surface 
observations and the results are used to compare with estimates based on radio- 
sonde observations or with estimates based on TOVS data. We find that the 
TOVS-based estimates of the stress are significantly more accurate than those 
found with a constant geostrophic drag coefficient, with a root-mean-square 
error about half as large. This improvement is due to stratification effects 
included in the boundary layer model. The errors in the sensible heat flux esti- 
mates, however, are large compared to the small mean values observed during 
the field experiment 
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1 . Introduction 

The surface fluxes of radiation, sensible and latent heat, and momentum are 
the links between the atmosphere and the ice. Over the polar ice pack there are 
few surface-based measurements of these quantities, and the measurements that 
do exist are for isolated points and are of relatively short duration or represent 
estimates from bulk measurements of wind and temperature (e.g., Leavitt, et al 
1978). Our understanding of the Arctic climate and our ability to predict the 
weather would be improved with accurate and timely estimates of these fluxes 
over the entire basin' Satellite observations hold some promise of determining 
these fluxes with sufficient accuracy. In this paper we explore the possibility of 
estimating surface fluxes of sensible heat and momentum using the TIROS 
Operational Vertical Sounder (TOVS) on the National Oceanic Atmospheric 
Administration (NOAA) series of polar orbiting satellites. Information retrieved 
from TOVS data has tremendous potential for studies of energy transfer and 
atmospheric dynamics over sea ice where conventional data arc sparse. This 
rich source has only begun to be tapped. Francis (1994). for example, uses 
retrievals to investigate the advection of heat from lower latitudes into and 
within the Arctic basin, and data from TOVS have also been used to estimate 
Arctic surface radiation fluxes (Schweiger and Key, in press). 

Large-scale models of sea ice require inputs of gcosirophic wind, down- 
welling short- and longwave fluxes and surface air temperature. Down welling 
short- and longwave fluxes are traditionally computed using simple parametcr- 
izations based on air temperature, cloud fraction, and humidity. Quesuonable 
climatologies or data from poorly validated numerical weather prediction mod- 
els are frequently used for these inputs. Tfrrbulent exchange of heat is specified 
using bulk formulae with a surface (skin) temperature determined from a sur- 
face energy balance equation. A constant geostrophic drag coefficient is used to 
determine the ice stress, ignoring changes in the boundary layer stratification 
that can have important ramifications for the actual stress (Overland and David- 
son, 1992). 

Our approach to the problem of estimating the turbulent fluxes using satel- 
lite data is to model temperature profiles in the ice and the atmosphere over 
time and to determine within the model the required sensible heat and momen- 
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turn fluxes needed to maintain the surface energy balance. We will use the 
observations from the LeadEx ice camp to verify the results. 

The numerical experiments described below demonstrate a simple one- 
dimensional thermodynamic sea ice model that is driven by input variables first 
derived from surface-based observations and then derived from satellites. The 
modeled atmospheric profiles of temperature and wind are allowed to vary in 
time so the surface stress can respond to the changing stratification of the 
boundary layer. The model is used to compute turbulent fluxes in a series of 
three experiments: 

Experiment 1, Surface Forcing: The first experiment is designed to estab- 
lish a baseline for comparison with the other two. It uses surface-based obser- 
vations only and does not include a boundary layer model. The surface 
turbulent fluxes calculated in this model are compared to the surface measure- 
ments. The forcing parameters are surface-observed air temperature, humidity, 
and winds for a height of 2 m, and downwelling short- and longwave radiative 
fluxes. The free parameters are the skin temperature, the ice temperature pro- 
file, and the surface turbulent fluxes. 

Experiment 2, Radiosonde Forcing: The second experiment is designed 
to show that the boundary layer portion of the model works correctly. The forc- 
ing parameters are radiosonde data (temperature, humidity, and wind vector) 
from the 850-mb level and surface-observed downwelling radiative fluxes. The 
free parameters are those in Experiment 1 plus boundary layer (850 mb to the 
surface) profiles of air temperature, humidity, and wind vector. 

Experiment 3, TOVS Forcing: The third experiment is designed to show 
the utility of satellite-based observations. The forcing parameters are the 
850-mb air temperature and humidity retrieved from TOVS, geostrophic winds 
from buoy-observed pressure fields (assumed to apply to the 850-mb level), and 
radiative flux estimates derived from TOVS radiances. Free parameters are the 
same as in Experiment 2. 

The forcings for the three experiments are summarized in Table 1 which 
includes the mean and standard deviation of each parameter. The fluxes com- 
puted in Experiments 2 and 3 are compared to those computed in Experiment 1. 
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2. Model Description 

We adapted a one-dimensional boundary layer model taken from the 
National Center for Atmospheric Research Community Climate Model 
Version 2 (NCAR-CCM2, Hack et al., 1993). We use only the components that 
compute the vertical diffusion, the surface energy balance, and the subsurface 
thermal flux. The model is essentially a column model that extends from the 
bottom of the ice to either the 2-m or 850-mb level in the atmosphere. The 
principal forcing parameters are air temperature and wind velocity at the forcing 
level (2 m or 850 mb) and surface downwelling radiative fluxes. Figure 1 
shows a schematic diagram of the model. 

We selected the NCAR-CCM2 model because it is well established in the 
scientific community, runs quickly, and can be forced with the observations we 
have at hand. Our adaptation of the CCM2 model leaves out the algorithm that 
determines radiative heating rates within the atmosphere. We will briefly dis- 
cuss each of the model components in turn, beginning with the ice model, then 
the surface energy balance model, and finally the boundary layer model. 

2.1 Thermodynamic ice model 

The thermodynamic ice model is in the tradition of Maykut and Unter- 
steiner (1971). We use a version of the CCM2 subsurface thermal model that 
includes enhanced ice physical properties, additional ice and snow layers, and 
penetration of solar flux below the surface. The temperature profile is calcu- 
lated by solving the thermal diffusivity equation for 10 layers with a fully 
implicit Crank-Nicholson scheme. Layer thickness, density, heat capacity, and 
thermal conductivity are specified. There are 5 quadraucally spaced layers in 
the snow and 5 in the ice. The top layer in the snow is just 1 cm thick, in order 
to resolve fluctuations on the order of one hour (Guest and Davidson. 1994). 
The solar radiative flux is allowed to penetrate the surface and be absorbed 
throughout the snow layer, with but a small fraction entering the ice (Grenfell. 
1979). Consistent with observations made at the LeadEx camp, the snow depth 
in the model is 10 cm, the snow density is 400 kg m'\ the ice thickness is 3.0 
m, and the albedo is 0.80. All four are held constant. 
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2.2 Surface energy balance 

The terms of the energy balance for the surface layer in the snow include 
the net radiation, F p the sensible heat flux, F s , the latent heat flux, F ql> and the 
conductive flux from below, F c . The balance is expressed as 

F t + F a + F q + F c = 0. (J) 

All fluxes are considered positive if directed toward the surface and the balance 
is zero because the top layer in the snow has negligible heat capacity and melt- 
ing does not occur. The net radiation consists of both the shortwave and long- 
wave component 

F r = (1 - a)F^ (1-/J+ eCF dlw - ar s!c 4 ), . (2) 

where F daw and F dlw are the downwelling shortwave and longwave fluxes at 
the surface, a is the albedo of the ice, I Q is the fraction of the net shortwave 
flux that penetrates below the first layer, e is the surface emissivity, a is the 
Boltzmann constant, and T s{c is the surface (skin) temperature. The down- 
welling radiative fluxes are prescribed either from direct measurements (Experi- 
ments 1 and 2) or from fluxes derived from TO VS radiances (Experiment 3). 

The sensible heat flux is determined by using a stability-dependent 
exchange coefficient (Hack et al., 1993), 

F s = pc p C s U(T tir -T s[c ), (3) 

where T tic is the model-calculated surface (skin) temperature, p is the air den- 
sity, and c p the heat capacity. The air temperature T& and the wind speed U is 
for a height of 2 m in Experiment 1 and the midpoint of the lowest layer (30 
m) in Experiments 2 and 3. The stress is also found using a bulk parameteriza- 
tion 

T = p C d t/ 2 . ' (4) 

The transfer coefficients C s and C d are calculated as functions of the bulk Rich- 
ardson number of the lowest layer. The latent heat is determined with the 
model-calculated surface (skin) mixing ratio and the mixing ratio determined 
either from the relative humidity measurements at 2 m (Experiment 1) or esti- 
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mated by the model for the lowest layer (Experiments 2 and 3). The transfer 
coefficient for humidity is the same as that for sensible heat. 

The 10-m neutral values of the heat transfer coefficient and the drag coeffi- 
cient are set to the average values measured at the camp, C d = C s = 0.0012. 
Note that this measured value of C g is corrected from that reported in Ruffieux 
et al. (1995). The uncertainty in the measured values of these coefficients is 
0.0002. The same coefficients are used in all three experiments. 

2.3 Atmospheric boundary layer model 

The boundary layer model is outlined in Holtslag and Boville (1993) and in 
Hack et al. (1993). It is an explicit, non-local parameterization of the vertical 
fluxes of heat, moisture, and momentum based on eddy diffusivities. The eddy 
diffusivities are the same for all three fluxes. The local diffusion scheme uses 
an eddy diffusivity determined independently at each point in the column, based 
on local vertical gradients of wind and virtual potential temperature, similar to 
the usual approach in global atmospheric models. The non-local scheme, used 
in unstable conditions only, determines an eddy diffusivity profile based on a 
diagnosed boundary layer height and a turbulent velocity scale. It can allow for 
a flux directed counter to the local gradient. The major advantage ot the non- 
local parameterization is that large eddy transport is accounted tor and entrain- 
ment effects are treated implicitly. 

There are 10 layers in the model. The first is centered at 'O m. the second 
at 75 m, and the rest are approximately 140 m thick. The boundary laser model 
provides tendency terms for temperature, moisture, and wind velocity at each 
level. Geostrophic adjustment terms are also included lor the wind velocity. 
The model time step is 15 minutes. 

3. Forcing parameters 

3. 1 Surface data 

The LeadEx ice camp occupied a muitiyear ice floe from 20 March to 22 
April 1992 in the Beaufort Sea about 400 km north of Deadhorse. Alaska. The 
purpose of the experiment was to study the physical processes associated with 
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leads in heavy pack ice. At the main camp a full set of meteorological observa- 
tions was made with 15-minute averaging periods. In addition, radiosonde 
ascents and measurements of radiative and turbulent fluxes were made. The 
surface observations, including the turbulent flux measurements, are reviewed 
by Ruffieux et al. (1995). The authors report significant diurnal cycles in 
observed skin temperature, air temperature, relative humidity, net radiation, and 
sensible heat flux. The average and standard deviation of each of the forcing 
parameters are shown in Table 1. Figure 2 shows the 2-m air temperature and 
the 2-m wind speed for the period 28 March (day 87) to 21 April (day 111), and 
Figure 3 shows the downwelling short- and longwave fluxes. The turbulent 
flux measurements for the period (sensible heat and stress) are shown in Figure 
4. Ruffieux et al. (1995) report that the measurement uncertainties are ± 2 
W nr 2 in the sensible heat flux, ± 0.003 Nt m -2 in the stress, and ± 3 W m -2 in 
the downward longwave radiative flux. The two shortwave flux instruments 
used in the field agreed to within 2% (about 8 W nr 2 at midday). 

3.2 Radiosondes 

The radiosonde observations from the LeadEx camp and nearby shore sta- 
tions are reviewed by Persson et al. (1992). They report that a well-mixed day- 
time boundary layer, with a depth of a few hundred meters, was common but 
not always present and that the median height of the inversion base was only 
28 m. The mediian inversion depth was 589 m. They also report that low-level 
jets were common and that the air in the ABL was often saturated with respect 
to ice. The ascents were made every 12 hours and hourly time senes for the 
850-mb level were generated with cubic spline interpolations. The air tempera- 
ture and wind speed for that level are shown in Figure 2. Note that the 850-mb 
temperature averages over 5°C warmer than the 2-m temperature. The mea- 
sured winds, both at the surface and aloft, indicate a high wind event in the last 
5 days of the period. This event is accompanied by significant warming as 
well. 

3.3 TOVS 

The TOVS instrument comprises three radiometer arrays, two of which are 
used in the retrieval. The High-resolution Infrared Radiation Sounder (HIRS) 
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has 24 channels in the visible through thermal infrared spectrum, and the 
Microwave Sounding Unit (MSU) has 4 channels in the 60 GHz oxygen band. 
Radiances were processed with a version of the Improved Initialization Inver- 
sion (31) algorithm (Chedin et al., 1985) that has been improved for applications 
in polar regions (Francis, 1994). Temperature profiles, humidity in 5 thick lay- 
ers, cloud information, and surface skin temperature are produced for each 
retrieval. Daily gridded fields of the various parameters at 100-km resolution 
north of 60°N are produced for the NASA Pathfinder Program and are available 
from the National Snow and Ice Data Center. This project uses individual 
retrievals, not the gridded fields. 

The accuracy of TOVS-derived air temperature usually decreases below 
900 mb, consequently we use the retrieved 850 mb air temperature to force the 
model. The wind speed at this level is assumed to be gcostrophic. determined 
from air pressure measurements made by surface buoys. There arc about 100 
TOVS retrievals each day within a 200-km range of the ice station. A smooth 
time series of the various parameters is derived from these retrievals by averag- 
ing simultaneously in time and space all observations within 200 km of the 
camp that were obtained within 1 day of the target time. The TOVS 850-mb air 
temperatures are compared to the radiosonde observations in Figure 2. The 
comparison is good except for the first 8 days when the TOVS temperature is 
up to 5 K warmer than the radiosonde measurements. This is a time of gener- 
ally greater cloud cover as shown in the higher values of the downwelling long- 
wave fluxes seen in Figure 3. The TOVS 850-mb air temperatures averaged 
1.5 K less than the radiosonde temperatures, with an RMS difference of 2.8 K 
and a correlation of R = 0.90. Also seen in Figure 2 is the ume senes of the 
geostrophic wind speed determined from air pressure measurements made by 
buoys and interpolated by the International Arctic Buoy Program. The mean 
geostrophic wind speed is 7.0 m s' 1 (Table 1). The 850-mb radiosonde winds 
average 0.6 m s' 1 more than the geostrophic winds and the RMS difference 
between them is 3.4 m s' 1 . 

Downwelling shortwave (DSW) and longwave (DLW) radiation fluxes at 
the surface are obtained directly from TOVS brightness temperatures using the 
method developed by Schweiger and Key (in press). This method (RADnet) is 
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based on a neural network that was trained using surface observations from the 
LeadEx experiment It uses TOVS HIRS and MSU brightness temperatures as 
inputs and provides estimates of DLW and DSW as outputs. 

The DLW fluxes were interpolated to hourly observations by averaging in 
time and space as discussed above. DSW observations, naturally displaying a 
large diurnal cycle, were interpolated to hourly values using a “solar zenith 
angle” weighting scheme. RADnet-computed DSW values at satellite overpass 
times were used to compute an “apparent cloudiness” by inverting the down- 
welling shortwave parameterization given by Shine (1984). This parameteriza- 
tion is designed to estimate the downwelling shortwave flux using the solar 
zenith angle, cloudiness, cloud optical depth, and surface albedo. For each 
retrieval, the cloud fraction required to obtain the TOVS-estimated DSW in the 
Shine parameterization is determined. The apparent cloudiness includes other 
radiative effects, such as DSW retrieval error, variable albedo, and variable 
cloud optical depth, as well as radiative effects not included in the parameter- 
ization. This apparent cloudiness is then averaged in time and space to form an 
hourly time series. The hourly cloud fractions are finally used in the Shine 
parameterization to compute hourly DSW fluxes. The mean differences 
between the measured downwelling fluxes and the TOVS-derived fluxes are less 
than 2 W nr 2 . The RMS differences are 17 W m~ 2 for DWL and 21 W m" 2 for 
DSW. The time series for both fluxes are compared to the surface-measured 
values in Figure 3 in which the DSW from TOVS is plotted only at local noon, 
since the TOVS-based estimate and the measured value from the camp are sim- 
ilar when viewed as time series. 

4. Results 

We use 24-hour running mean values for numerical comparisons applied to 
both the measured and the modeled fluxes. We do this because the forcing val- 
ues for both the radiosondes and for TOVS are derived from either twice daily 
samples (radiosondes) or from running 24-hour mean values (TOVS). In addi- 
tion, variations on scales of less than 1 day are of limited interest for most cli- 
mate studies. The mean and standard deviations of the measured turbulent 
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fluxes are shown in Table 2, in which the standard deviation is indicated for 
hourly averages and daily averages. Times when the wind direction corre- 
sponded to the direction of the camp have been excluded. Comparisons of the 
surface-flux values from the three numerical experiments are summarized in 
Table 3. First, Experiment 1 is compared to the measurements, then Experi- 
ments 2 and 3 are compared to results from Experiment 1 . 

Figure 4 shows time series of sensible heat flux and stress for Experiment 1 
compared to hourly averages of the surface measurements. The diurnal cycle of 
the sensible heat flux is well represented, but significant biases appear in the 
last few days when the winds are strongest. The change in the magnitude of 
the differences may be associated with unknown changes in the snow depth, 
density, or albedo that are associated with the wind storm or with changes in the 
measurement systems. Table 3 compares the daily (24-hour running mean) val- 
ues of the results of Experiment 1 to the measurements. The bias is -I.3W m' 2 . 
The RMS difference is larger, 6.2 W m -2 , reflecting errors in the model, includ- 
ing characterization of the surface, or errors in the flux measurements. Errors in 
either the turbulent flux measurements or the radiative flux measurements would 
both contribute to differences between the measured and the modeled fluxes. 
The average model stress is 0.011 Nt nT 2 less than the measured stress. A 
value of C d = 0.014 provides near perfect agreement, a value that is within the 
range of uncertainty in the reported value (Ruffieux et al. 1995). 

Experiments 2 and 3, which include the boundary la\er. are compared to 
the results of Experiment 1 in Figure 5 and in Table 3. In the radiosonde- 
based estimates of the sensible heat flux. Experiment 2. the bias and the RMS 
difference are 3.8 W m -2 and 11.8 W nr 2 , respectively. With the TOVS-based 
estimates, Experiment 3, these differences drop slightly to 3.0 W in'- and 
9.1 W m -2 . The stress estimates are also slightly better with the TOVS-based 
model. The bias is 0.016 Nt nO 2 for Experiment 2 and -0.003 Nt in " lor 
Experiment 3. The RMS differences arc 0.052 Nt m - and 0.021 Nt m ~ 
respectively. The improved performance of the TOVS-based model may be due 
to the increased averaging (horizontal, vertical, and temporal) inherent in the 
TOVS-based forcing. It may also reflect the presence of significant thermal 
wind effects which are not included in the model. These effects would make 
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the 850-mb radiosonde wind a poorer basis for estimating the surface stress 
than the surface geostrophic wind that is used in Experiment 3. 

5. Discussion and Conclusions 

One of our principal motivations for this study is to test a scheme in which 
estimates of the surface stress can account for variations in the boundary layer 
stratification. How satisfactory are traditional approaches, in which stratifica- 
tion is not accounted for? The stress in large-scale sea-ice models is often cal- 
culated with just the geostrophic wind speed G and a fixed geostrophic drag 
coefficient C g as 

X = pC l G 1 (5) 

In our numerical experiments we have used a measured value of the local 10-m 
drag coefficient and so, to provide a useful comparison, we select a value of C g 
that gives zero bias in the stress calculated from (5) when compared to the 
Experiment 1 results. That value is C g = 0.00051. The RMS difference 
between the Experiment 1 estimates and the estimate from (5) is 0.033 Nt nr 2 
compared to 0.021 Nt m' 2 for the TOVS-based estimate (Experiment 3) which 
also showed negligible bias. Thus the error in the TOVS-based estimate is 
36% less than that found with a constant geostrophic drag coefficient. This 
improvement is due to the inclusion of stability effects in the boundary layer 
model that influence the geostrophic drag coefficient and are determined by the 
surface energy balance. 

In these numerical experiments, the satellite-based estimates of the stress 
are significantly better than estimates using a constant drag coefficient. Because 
of this, we feel that a column model with TOVS-based forcing shows promise 
for more accurately determining the surface stress. The heat- flux estimates, on 
the other hand, show errors of several watts per meter squared. Only the high 
wind event near the end of the period is detected, and even that is with poor 
accuracy. The heat- flux estimates from TOVS-based forcing (Table 3) are, 
however, as accurate as those from radiosondes. 

Experiment 1 highlights some of the difficulties in using an energy balance 
model to determine the sensible heat flux. First, the net radiation must be deter- 
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mined, which depends both on accurate downwelling fluxes and on accurate 
estimates of the albedo, and, second, the conductive flux must be determined, 
which depends on an accurate estimate of the snow depth and snow density. 
Finally, the fluxes during LeadEx are small, often averaging over the course of 
a day to less than 10 W nr 2 in absolute value, so measuring a significant por- 
tion of their variance is difficult. Space-based estimates in the near future will 
probably be able to detect only the strong events such as that seen in the last 6 
days. Although this event is detected in the TOVS-based estimates (Figure 5), 
significant bias remains. 

The results are also sensitive to the turbulent transfer coefficients specified. 
We found that if the drag coefficient and the heat transfer coefficient were tuned 
independently to provide the minimum error in Experiment 1. the RMS differ- 
ences between the results of Experiments 2 and 3 and the measurements could 
be made to fall below 5 W m -2 and the bias in the stress compared to the mea- 
surements made to disappear. An accurate characterization of the surface 
roughness will be needed to adequately determine the stress. 

These experiments have shown that it is feasible to estimate the surtaec tur- 
bulent fluxes using data collected from satellites. However much work remains 
to be done before our procedure can be coupled to a two-dimensional ice/ocean 
model and fluxes can be computed on a routine basis for the entire Arctic. 
First, additional work needs to be done to select and tune the appropriate 
numerical model. Is the NCAR CCM2 the right one? Other models may prove 
more useful. Second, the results of the model are sensitive to a number ot dif- 
ferent parameters including the albedo, snow depth, ice thickness, and the tur- 
bulent exchange coefficients. Determining the correct values to use in the 
absence of in situ data could be a formidable task. Finally, additional valida- 
tion experiments should be performed for different seasons of the year and dif- 
ferent ice conditions. 

Possible enhancements to the model might include the addition ot a radia- 
tive flux divergence estimate and lead heat flux estimates. Both would contrib- 
ute to the heating rates within the boundary layer. The inclusion of clouds 
within the boundary layer would improve the flux divergence calculations. The 


IGS ‘96. TOYS Surface Fluxes 


September 4, 1996 


12 


sensible heat flux determined here is appropriate for thick ice. The method 
could be extended to other ice classes and to estimates of regional heat flux by 
estimating the near-surface air temperature and wind speed for the dominant ice 
class in a region using the TOVS-based forcing and the column model 
described here. The near-surface air temperature and wind speed could then be 
used to estimate the heat fluxes from other ice classes, including leads. 
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Table 1: Forcings for Three Numerical Experiments 


Experiment 

Parameter 

Mean 

Standard 

Deviation 

1 : Surface Forcing 

2-m Air Temperature (°C) 

-22.8 

4.2 


2-m Wind Speed (m s" 1 ) 

3.8 

2.6 


Downwelling Longwave (W m -2 ) 

180.3 

27.4 


Downwelling Shortwave (W m~ 2 ) 

140.4 

150.6 

2: Radiosonde Forcing 

850-mb Air Temperature (°C) 

-16.9 

5.4 


850-mb Wind Speed (m s" 1 ) 

7.6 

5.9 


Surface-measured Downwelling 
Longwave (W m -2 ) 

180.3 

27.4 


Surface-measured Downwelling 
Shortwave (W m -2 ) 

140:4 

150.6 

3: TOYS Forcing 

850-mb Air Temperature (C) 

-18.4 

■ 4.3 


Geostrophic Wind Speed (ms' 1 ) 

7.0 

5.1 


RADnet Downwelling 
Longwave (W m~ 2 ) 

178.7 

19.9 


RADnet Downwelling 
Shortwave (W nr 2 ) 

142.0 

150.0 


Table 2: Measured Thrbulent Fluxes 


Flux 

Mean 

Standard Deviation, 
hourly averages 

Standard Deviation, 
daily averages 

Sensible Heat (W m' 2 ) 

2.9 

10.5 

7.2 

Stress (Nt m -2 ) 

0.070 

0.090 

0.083 


IGS '96. TOYS Surface Fluxes 


September 4, 1996 


15 


Table 3: Sensible Heat Flux and Stress Comparisons (Daily Averages) 


Experiment 

Number of 
Days 

Bias 

RMS 

Difference 

Difference 

Standard 

Deviation 

1: Surface forcing, 

compared to measurements 





Sensible Heat (W m -2 ) 

17 

-1.8 

6.2 

5.9 

Stress (Nt m -2 ) 

17 

-0.011 

0.021 

0.018 

2: Radiosonde forcing, 
compared to Experiment 1 





Sensible Heat (W m~ 2 ) 

24 

3.8 

11.8 

11.1 

Stress (Nt m -2 ) 

24 

0.016 

0.052 

0.050 

3: TOVS-based forcing, 
compared to Experiment 1 





Sensible Heat (W m -2 ) 

24 

3.0 

9.1 

8.6 

Stress (Nt m -2 ) 

24 

- 0.003 

0.021 

0.021 

Geostrophic wind only, 
compared to Experiment 1 





Stress (Nt m -2 ) 

24 

0.000 

0.033 

0.033 


rfr 
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List of Figures 


Fig. 1. Schematic of the model temperature profile and heat 
fluxes. The forcing parameters are the surface downwelling 
radiative fluxes and the air temperature and wind speed at either a 
height of 2 m (Experiment 1) or at the 850-mb level (Experiments 
2 and 3). There are 5 layers within the snow, 5 within the ice, and 
10 within the atmospheric boundary layer. 

Fig. 2. Comparisons of the forcing air temperatures (a) and the 
forcing wind speeds (b). 

Fig. 3. Comparison of the measured downwelling radiative 
fluxes at the LeadEx camp to those derived from the TO VS 
radiances: (a) shortwave and (b) longwave. The shortwave fluxes 
are nearly identical because of the common solar-zenith-angle 
dependence, so only the daily maxima in the TOVS-based 
estimates are plotted. 

Fig. 4. Surface turbulent fluxes from Experiment 1 with surface 
forcing (line) and from measurements (points): (a) sensible heat 
flux, (b) stress. The measurements are hourly averages and the 
model output is hourly. 

Fig. 5. Comparisons of the three numerical experiments: (a) 
sensible heat flux, (b) stress. The time series are 24-hour running 
means of the model outputs. 
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Figure 1. Schematic of the model temperature profile and heat fluxes. 
The forcing parameters are the surface downwelling radiative fluxes and 
the air temperature and wind speed at either a height of 2 m (Experiment 1) 
or at the 850-mb level (Experiments 2 and 3). There are 5 layers within the 
snow, 5 within the ice, and 10 within the atmospheric boundary layer. 
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Exp. 1: Surface Forcing 
Exp. 2: Radiosonde Forcing 
Exp. 3: TOVS Forcing 


Exp. 1: Surface Forcing 
Exp. 2: Radiosonde Forcing 
Exp. 3: TOVS Forcing 


Figure 5. Comparisons of the three numerical experiments: (a) sensible 
heat flux, (b) stress. The time series are 24-hour running means of the 
model outputs. 
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